Gravitational arcs as a perturbation of the perfect ring. 
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ABSTRACT 

The image of a point situated at the center of a circularly symmetric potential 
is a perfect circle. The perturbative effect of non-symmetrical potential terms 
is to displace and break the perfect circle. These 2 effects, displacement and 
breaking are directly related to the Taylor expansion of the perturbation at first 
order on the circle. The numerical accuracy of this perturbative approach is 
tested in the case of an elliptical potential with a core radius. The contour 
of the images and the caustics lines are well re-produced by the perturbative 
approach. These results suggests that the modeling of arcs, and in particular of 
tangential arcs may be simplified by using a general perturbative representation 
on the circle. An interesting feature of the perturbative approach, is that the 
equation of the caustic line depends only on the values on the circle of the lens 
displacement field along the 9 direction. 

Subject headings: gravitational lensing:strong lensing 

1. Introduction. 

Since the discovery of gravitational arcs in clusters of galaxies by Linds & Petrosian 
(1986) and subsequently by Soucail et al. (1987), the observation and study of arcs has 
developed considerably, and is now becoming an essential tool in astrophysics. Arcs provides 
a wealth of information about the mass distribution in clusters of galaxies (see for instance: 
Comerford et al. 2006, & Broadhurst el al. 2006). However, the mass distribution of the 
astrophysical lenses is complex and involves a large parameter space which is difficult to 
explore. Thus, the derivation of a simplified perturbative theory that is able to re-produce 
the general features of gravitational arcs is an interesting tool to help understanding complex 
gravitational lenses. 
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2. Basic ideas 

Let's assume that the projected density of the lens ps is circularly symmetric and cen- 
tered at the origin. Let's also assume that the lens is dense enough to reach critical density 
at a given radius Re- Under such hypothesis, the image by the lens of a point source placed 
at the origin will be a perfect ring. We are now interested in small perturbations of this 
perfect ring. There are two type of perturbations to the point source perfectly aligned in a 
circularly symmetric potential: first the source may not be perfectly at the center, and sec- 
ond the potential may not be perfectly circular. Let's be more specific, in polar coordinates, 
the lens equation writes: 

rs = r - Ur - -TTT ue (1) 



dr J \r dd 
In the unperturbed case, the equation reads: 

r-#^=0 (2) 
dr 

Where 0o is a function of r only. Let's now perturb this equation by introducing a small 
displacement of the source from the origin: Tg, and a non-circular perturbation to the po- 
tential, ^. Note that the perturbation on and the perturbation on cj) are assumed to be 
of the same order. The perturbation may be described by the following formula: 

(3) 

Here e is a small number: e << 1. Given a position (r^) for the source, the image posi- 
tions (r, ^), can be obtained by solving equation ([T]). However solving directly may prove 
impossible in the general case. It is easier to find a perturbative solution by inserting Eq. 
([3]) in Eq. ([1]). Assuming that e is small, the perturbation will introduce a deviation from 
the perfect circle that will be of order e. An interesting feature of the un-perturbed solution 
is that the image of a single point at the origin is a full circle, which covers all range in 
9. Thus whatever the position 9 of the perturbed solution, there is always a point at the 
same 9 in the un-perturbed solution. However, the point in the un-perturbed solution will 
be located at a slightly different radius r. The response to the perturbation on r may be 
written, r = Re + e dr. For convenience, it is always possible to re-scale the coordinate 
system, so that the Einstein radius is exactly equal to unity, in this case, the perturbation 
on r is: 

r = 1 + e dr (4) 

To summarize, we must solve perturbatively Eq. ([1]), by expanding around the un-perturbed 
solution (1, 6*) for small values of e. Note that this requires to expand the potential at r = 1. 
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Using Eq. ([3]), the Taylor expansion of (j) may be written: 

oo 

(t) = ct>, + ei, = Y,[Cn + e Um (r - 1) 



Where we define, 



And: 



Cn 



n=0 
1 



(r=l) 



(r=l) 



(5) 



(6) 



(7) 



It is now possible to expand each side of Eq. ([T]) in series of e, in the vicinity of the un- 
perturbed solution. By inserting Eq. ([5]), and Eq. (jl]) in Eq. ([1]), the response dr to the 
perturbation defined in Eq. ([3]) can be estimated to the first order in e. 

e Ts = {I + e dr - Ci - e 2 C2 dr - e fi) Ur - e — - 

ot) 



Note that Eq. ([2]) implies that Ci = 1, consequently: 



Fs = {k dr - fi) Ur - — Ug 



With: K = 1 - 2 Co 



3. Reconstruction of images. 

3.1. Circular source contours. 

Let's consider a circular contour on a source with center ro radius, Rq. The equation 
for this contour is: 

(r, - ro)' = Rl (9) 

Note that effect of the translation by the vector ro = {xQ,yo) can be taken into account by 
re-defining /o, and fi in Eq. (IT^ : 



rs = [k dr- /i) Ur 



dfo 



_fi = fi + xo cos 9 + yo sine i = 0, 1 
The image by the lens of this contour can be obtained using Eq's. and (ITUil : 



(10) 



(11) 
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Solving Eq. ffTTl) . the 2 following solutions for dr can be obtained: 



1 




(12) 



In the above equation the eondition for image formation i. A = fl? - (f) > 0. The 

mean position of the 2 contour lines is — , and the image width along the radial direction is 

2 A 



3.2. Elliptical source contours. 



Using Eq. ffTOj) the source Cartesian coordinates writes: 

^= (^ndr- /i) COS0 + ^sin^ 
Us = [k dr - /i) sin 6* - ^ cos 6 



(13) 



The equation for an elliptical contour aligned with its main axis aligned with the coordinate 
system is: 

{l-r^)xl + {l + r])y'^=Rl (14) 
Using Eq.'s (IT^ and (IT^ one obtains: 




rj sin 26 dfo i 



^Rl S-[1W) (^)' 



(15) 



3.3. Numerical testing 



The accuracy of Eq. f|T2l) has been tested for a source with a circular contour by 
direct comparison with numerical integration of the lens equation using ray-tracing. For this 
application, the potential proposed by Blandford & Kochanek (1987) (hereafter, BK1987) 
was chosen. The asymptotic isothermal form of the potential was selected. The potential 
reads: 



s vs^ 



+ 1 



1 + 



(1 -r/cos20) 



- 1 



(16) 



The ellipticity of the potential is controlled by the parameter ?], BK1987 conclude that for 
physical reasons rj must be less than 0.2. The circular, unperturbed potential 0o corresponds 
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to r/ = 0, in this case, at r = 1, = 1 which imphes that Re 
and (I7j) the parameters of the perturbative expansion are: 



1. According to Eq's. 



' d/b _ r? sin 26* 

-y/s2+l_r;cos2e 



/l 



Vs^ + l{l-T] COS 20) _ -j^ 

y/s'^+l-rjcos 26 



(17) 



Using equations (IT^ combined with Eq. (fT7|) it is possible to trace the 2 hnes that describe 
the image contour. The 2 solutions for dr join at the contour edges. The result is visible in 
Fig. (1). Note that the large tangential arc is well reproduced by the approximation while 
there is some offset in the position of the inner image. The problem of the inner image 
can be corrected using an iterative approach: the perturbative method gives a first guess 
of the image position, at this initial position one may carry a local Taylor expansion of the 
potential in order to find a better solution, and the procedure may be iterated again. 



3.4. Inverse modelling 

Given a system of arcs or arclets it is possible to re-construct the perturbation field. 
Let's define the arc system to reconstruct as a set of contours, with one contour per image. 
For each contour a radial line of direction 6 intersect the contour in 2 points, ri and r2 
(provided that 6* is in a suitable range). It is simple to relate these 2 functionals to the fields, 
/i and /o of the perturbative approach, in the case of a source with circular contour (Eq. 

my- 



18) 



Where C is a constant term. Note that k is unknown, actually there is a fundamental 
degeneracy in k, only fi/n, fo/i^, and Rq/k may be determined. By using inner of radial 
image in the iterative approach mentioned in Sec. (3.3), in some case it might be possible to 
break the degeneracy on k. It is clear also that instead of circular contours, one could have 
considered elliptical one, and Eq. (ITSIl may have been used. In such case, it is just a matter 
of where to put the complexity, either in the source or the lens, and obviously the judging 
criteria should be that as a whole the complexity is minimal. 
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4. Caustics in the perturbative approach. 

Caustics are singularities, which are defined by the simple property that the determinant 
of the Jacobian Matrix J is zero on the caustic lines. 

J _ dxs_dy^ _ dxs_dys_ _ ^ 
dr do do dr 

The calculation of the Jacobian is straightforward from Eq. (fT3|) . it follows that: 



dr = — 

K 



d^O 



(19) 



Note that in this case a shift by a vector tq does not have to intorduced thus, fo = fo, and 
/i = /i. Eq. fin?]) defines the critical lines. The caustics in the source plane can be obtained 
by inserting Eq. (fT9|) in Eq. (fT3|) : 

x, = ^cos^ + ^sin^ 
^ = 4§^sin ^-f cos^ 

Not that Eq. ( l20i) line depends only on fo, which is directly related to the multipole expansion 
of the potential on the circle. The multipole expansion has also the advantage to relate 
directly /o to the density by the means of the coefficients of the multipole expansion at 
r = 1. Turing now to a numerical application, by inserting the expression of /o given in Eq. 
( ITTll in Eq. (l20l) a parametric equation of the caustic line is obtained, the result is presented 
in Fig. (2). The perturbative calculation of the caustic curve is accurate, this suggest that 
this approximation may be used to derive general results on caustics. A simple result is that: 

r'=(^X + ^ (21) 



^dO J d^O 

Eq. (12T]) may be integrated other 0, the second term due to periodicity, it follows that the 
typical size of the caustics is directly related to the variance of This result demonstrates 
that the caustics cross-section, are closely related to the deviations of the potential to circular 
symmetry, which is confirmed by the numerical analysis of Meneghetti et al. (2007). 



The author would like to thank J. P. Beaulieu and S. Colombi for reading this paper. 
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Fig. 1. — A circular source of radius Rq = v^/20 Re is placed at [xq = ^,yo = 0] from 
the potential center along the X axis. The core radius of the potential is s = and 
the ellipticity is 77 = ^. The dark areas corresponds to the numerical solution obtained 
by ray tracing. The outer contour of the source is reconstructed using Eq's f|T2|) and f|T7j) 
and super-imposed (as a red line) on the ray tracing solution. The red contour is close to 
the outer contour of the ray-tracing solution for the tangential image. Due to the strong 
non-linearity of the potential near the center, the approximation is not as good for the inner 
image, although as mentioned in Sec. (3.3) this may be corrected using an iterative approach. 




Fig. 2. — Caustics of the BK1987 potential, the perturbative solution is plotted in black 
next to the numerical solution of the system of equations without approximations. Note the 
closeness of the 2 solutions. 
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